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We present a method for recovering the distribution functions of edge-on thin 
axisymmetric disks directly from their observable kinematic properties. The most 
generally observable properties of such a stellar system are the line-of-sight velocity 
distributions of the stars at different projected radii along the galaxy. If the gravita- 
tional potential is known, then the general two-integral distribution function can be 
reconstructed using the shapes of the high-velocity tails of these line-of-sight distribu- 
tions. If the wrong gravitational potential is adopted, then a distribution function can 
still be constructed using this technique, but the low-velocity parts of the observed 
velocity distributions will not be reproduced by the derived dynamical model. Thus, 
the gravitational potential is also tightly constrained by the observed kinematics. 

Key words: Galaxies: structure - Galaxies: kinematics and dynamics 



1 INTRODUCTION 

The dynamics of a collisionless stellar system is fully spec- 
ified by its distribution function (DF) - the phase density 
of its constituent stars - and the form of the gravitational 
potential that binds it together. The most fundamental goal 
of galactic dynamics is the recovery of this information from 
the observable properties of a galaxy. There is some hope 
that this objective might soon be achieved: sophisticated 
data analysis techniques applied to the Doppler broadening 
of galaxy spectra now routinely yield reasonable estimates 
for the full distribution of the stars' line-of-sight velocities, 
generally referred to as their line-of-sight velocity distribu- 
tion, or LOSVD (see e.g. Franx & Illingworth 1988, Bender 
1990, Rix & White 1992, van der Marel & Franx 1993, Win- 
sall & Freeman 1993, Kuijken & Merrifield 1993). In princi- 
ple, one can therefore measure a complete three-dimensional 
function, the density of stars as a function of both their line- 
of-sight velocities and the two spatial coordinates giving the 
position on the plane of the sky. Since the DF of a stellar 
system in equilibrium depends on at most three isolating 
integrals of motion (Binney & Tremaine 1987), it is plausi- 
ble that one might deproject the observed three-dimensional 
function into the galaxy's intrinsic DF. 

Although this dimensional argument suggests that the in- 
version might be carried out, we have no guarantee that the 
solution will be unique: it is possible that more than one 
DF may produce identical LOSVDs. Further, the dynamics 
of a galaxy depends on the form of the gravitational po- 
tential that confines it. It is therefore also possible that a 
given set of kinematic observations will be consistent with 
different dynamical models depending on the form adopted 



for the gravitational potential. Thus, we do not yet know 
whether the dynamics of a galaxy is uniquely specified by 
its observable properties. 

In some special cases, we can make more definitive state- 
ments. For a spherical system with a known gravitational 
potential, for example, it has been proved that the dynam- 
ics is completely specified by the observable kinematics. De- 
jonghe & Merritt (1992) considered the moments of the ve- 
locity distribution - the velocity dispersion and its higher 
order analogues. For each order, they demonstrated that 
the intrinsic moments of a spherical galaxy's velocity dis- 
tribution could be inferred from the observable line-of-sight 
component. They thus elegantly proved that all the intrin- 
sic dynamical properties of the galaxy are specified by ob- 
servable quantities. Dejonghe & Merritt also showed that 
the observable kinematics will limit the possible forms of 
the gravitational potentials that might be confining the sys- 
tem, although not necessarily to the point of specifying it 
uniquely. 

A further significant advance was made by Merritt (1996), 
who considered edge-on axisymmetric systems where the DF 
respects two integrals of motion. He showed that for such 
systems one can use the information in the projected den- 
sity of stars and the first two moments of their line-of-sight 
velocity distributions to estimate both the DF and the grav- 
itational potential. Although not constituting a formal proof 
of uniqueness, the numerical solutions to Merritt 's equations 
do seem to indicate that accurate dynamical models can be 
inferred from the observable kinematic properties. 

In theoretical studies of galaxy dynamics, manipulation of 
the moments of the velocity distribution is mathematically 
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attractive because these quantities obey the Jeans equations 
and their higher-order analogues (Merrifield & Kent 1990, 
Dejonghe & Merritt 1992) . However, the application of such 
analyses to observations is, in practice, fraught with difficul- 
ties. The higher-order moments on which Dejonghe & Mer- 
ritt's (1992) proof depend are extremely sensitive to noise, so 
cannot be calculated from real data. Even the velocity dis- 
persion, calculated from the second moment of the LOSVD, 
is difficult to estimate from noisy data, particularly if the 
velocity distribution has long tails out to high velocities. 
In poorer-quality data, the velocity dispersion is often es- 
timated by fitting a Gaussian to the velocity distribution; 
there is in reality little justification for such a procedure, 
and systematically-biassed estimates for the dispersion are 
almost certain to be the result. 

In this paper, we therefore adopt a somewhat different ap- 
proach. The prevalence of techniques for measuring full 
LOSVDs suggests that it might be worth moving away from 
analysis of moments to using the complete shape of the ve- 
locity distribution directly in the analysis. To this end, we 
have been investigating the properties of a thin axisymme- 
tric disk model. In Section 2, we discuss how the dynamics 
of an inclined thin disk galaxy can be inferred directly from 
its observable kinematics. Such a simplified model is clearly 
a rather special case, since each line of sight passes through 
the disk at only a single point. Therefore, in Section 3 we 
consider the extra complexity that is introduced when such 
a disk is viewed edge-on, so that each line of sight passes 
through many points in the galaxy. We have used the ana- 
lysis of the inclined case to provide the basis for an iterative 
scheme for recovering the DF of such an edge-on disk. In 
Section 4, we test this algorithm on a realistic disk model, 
and show how the DF can, indeed, be recovered from the 
observable kinematics. We further show that even small er- 
rors in the assumed gravitational potential lead to inconsis- 
tencies in the inferred kinematic properties, and hence that 
the potential is very tightly constrained by the observable 
properties of such a galaxy. Section 5 discusses the possi- 
ble practical application of this analysis to studies of disk 
galaxies. 



2 KINEMATICS OF THIN DISKS 

The simplest plausible model for a disk galaxy is one in 
which the system is axisymmetric, with all the motion con- 
fined to a single plane. For such a system, Jeans' theorem 
states that, in equilibrium, the DF can depend only on two 
isolating integrals of motion, the energy of stars in the disk 
plane, and the angular momentum of stars about the centre 
of the galaxy (Binney & Tremaine 1987) . Defining the usual 
polar coordinates about the centre of the system and wri- 
ting the gravitational potential in the plane of the galaxy as 
$?(r), these integrals can be written as E = ~(«J?+«2)+\E , (r), 
and L = rw^, and the DF can be expressed as a function 
f{E,L). 

At any radius in the galaxy, the tangential velocity distribu- 
tion can be obtained by integrating the DF over the whole 
range of radial velocity: 



where 



f(E, L) dv r = 2 



f(E,L) dE 



If we transform variables and write 
U(r,v<f) = f4>(W<p,L), 



(2) 



(3) 



then equation (|l|) is an Abel integral equation relating 
f(E, L) and f$(W$, L). Thus, if f<j> is known, one can use a 
standard Abel inversion to solve for the full DF, 



f(E,L) = -~ 
n 



dU{W^,L) dW^ 

dW* y/2(W^-E)' 



(4) 



For a thin disk at inclination i to the line of sight, a quan- 
tity related to f$ that one can observe is the line-of-sight 
velocity distribution as a function of projected radius along 
the galaxy's apparent major axis, F ma j(r p ,vi os ). Via some 
simple geometry, one can show that 



Ejnaj y^p , ^los ) 



(r = r p ,v$ = v los / sin i)/ sin z. 



(5) 



Thus, as first pointed out by Merrifield & Kuijken (1994), 
if one observes the LOSVD along the major axis of such a 
disk, then one can readily derive /</, and hence solve for the 
DF by inverting equation (|l|) . More recently, a sophisticated 
investigation of this inversion has been made by Pichon & 
Thiebaut (1998), who showed that it is possible to regularize 
a non-parametric algorithm to carry out the requisite Abel 
inversion without excessive noise amplification. 
There is a degree of redundancy in the inversion, which 
means that its implementation can be simplified somewhat 
further. Every pair of {W$,L} values corresponds to two 
pairs of {r,v$} values, and hence two points in the obser- 
vable {r p ,vi„ s } parameter space. Thus, fy(W<t„L) is over- 
specified by the observed F ma j(r p ,vi oa ). If we define the 
galaxy's circular speed, v c (r), by the usual relation, 

%{r)=r—, (6) 

then it is straightforward to show that the DF is com- 
pletely specified by both the part of F ma j(r p ,Vi os ) where 
|«; s| > v c (r p ) sin i and by the part where \vi os \ < v c (r p )smi 
(Merrifield & Kuijken 1994). 

Using just one part of the LOSVD allows us to simplify 
the quadrature in equation (Q) . For the high- velocity side of 
the LOSVD in a realistic disk, f${W<p, L) is a monotonically 
decreasing function of ]¥</, at a given value of the angular 
momentum L. One can show that any positive distribution 
function such that 



df(E,L) 
dE 



< 



(7) 



(1) 



yields a tangential velocity distribution f$(W$,L) with the 
above property. In stellar dynamics one usually considers 
systems with distribution functions that obey equation (R). 
Definite criteria for the stability of spherical and axisym- 
metric stellar systems have been established for distribu- 
tion functions which are decreasing functions of the energy 
(Antonov 1962). Moreover it is found that most distribution 
functions used to produce realistic models of galaxies do ac- 
tually obey that condition. On the other hand, very little 
work, either numerical or theoretical, has been done using 
distribution functions which do not satisfy such a condition 
and very little is known about the stability of these systems 
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(Henon 1973, Perez & Aly 1996). Therefore, in our paper, 
we follow the generally accepted assumption given by equa- 
tion (Q). In this case can be written as a monotonic 
function of f^>(W^, L) at given L. This property allows us to 
transform the above relation into 



f(E,L) 



U(.W4=E,L) 



y/2(W4f'^,L))-E) 



(8) 



Provided that the function fj,(W<f,,L) can be numerically 
inverted to give W^(f^, L), this equation means that the DF 
can be calculated with a simple one-dimensional integration, 
bypassing the computationally-unstable estimation of the 
derivative of f^(W^,L) in equation (0). 



3 KINEMATICS OF EDGE-ON THIN DISKS 

As mentioned in the Introduction, an inclined thin disk is 
rather simpler than most realistic galaxy models, since any 
line of sight will only intersect it at a single point. Further, 
for an inclined system, one cannot tell from the photometry 
whether the assumption that the disk is thin is a valid one - 
it is only when the system is viewed edge-on that its razor- 
thinness will be unambiguously apparent. We therefore now 
turn to consider the case of such a system viewed edge-on. 
In this case, the LOSVD is not just a rescaled version of 
/(/.(W 7 ^, L) as in equation but instead is generated by an 
integration along the line of sight through the entire disk. By 
choosing the z-axis as the direction of the line of sight, one 
can write the LOSVD at some distance r p from the centre 



of the galaxy as 



F(r p ,Vi os ) 



dz \ dv±f(E,L) 



(9) 



where (vi os ,v±) are the components of the velocity parallel 
and perpendicular to the direction of the line of sight re- 
spectively. The two components of the velocity (v r ,Vj,) can 
be written in terms of (vi os ,v±) as: 



v r = v± cos <p + vi oa sm 4 
= vios cos (f) -v ± sin 4 

with 

cos <p = — ; sm <j> = — ; 
r r 

The integrals of motion are then 

e = i(«i + «L) + *(0 

L = rv<p = r p vi os — zv±. 



= y/rp + z 2 



(10) 
(11) 

(12) 

(13) 
(14) 



The extra integration along the line of sight clearly compli- 
cates the relationship between the observed kinematics and 
the intrinsic dynamics of the galaxy. However, the properties 
that we have derived above for inclined disks suggest a way 
to proceed. In particular, the fact that the DF is completely 
specified by the high- velocity tail of f$ can be used to our ad- 
vantage. In general, the mean-streaming circular motions of 
the stars in a realistic disk will be significantly greater than 
their random motions, so one can think of stars at each ra- 
dius moving with a mean streaming velocity of close to the 
local circular speed, with a smaller amount of random mo- 
tion superimposed. Thus, the LOSVD at a projected radius 



r p in an edge-on disk will be peaked at a velocity close to 
v c (r p ). At lower line-of-sight velocities, there will be contri- 
butions both from stars that have intrinsically low velocities 
and from stars that lie at large radii in the galaxy, so that 
their large circular motions are oriented mostly transverse 
to the line of sight, resulting in a small line-of-sight compo- 
nent. However, at line-of-sight velocities greater than v c (r p ), 
the majority of the stars contributing to the LOSVD will be 
those where the circular motion is oriented along the line- 
of-sight, which occurs only for stars at radii r ~ r p . Further, 
for the stars at radii close to r p , the line-of-sight component 
of their random motions measures their velocities in the 4> 
direction. Thus, to a simple approximation, for vi os > v c , 



F(r p ,vi oa 



v toa ) x A, 



(15) 



where A is the length of path through the galaxy that lies 
sufficiently close to r = r p for there to be a significant con- 
tribution to F(r p ,vi os ). In terms used in the study of the 
Milky Way, this path-length covers the region of the tan- 
gent point, over which f<j> will not change rapidly (since r is 
not changing significantly), so the integral along the line of 
sight can be replaced by the simple product of equation ( |l5| ) . 
Thus, equation ( |l5| ) implies that by observing the part of 
the LOSVD at high velocities, we obtain an estimate for 
the shape of the high-velocity part of the tangential velocity 
distribution, which we need to estimate the DF. 
The close relation between F(r p , vi os ) and /^(r, v$) provides 
a good indication that the inversion of equation may 
not be significantly more ill-conditioned than was the case 
for the simpler case of the inclined disk [equation (Q)]. Al- 
though equation (Ji^) contains an extra integration, which 
often serves to smooth out some of the information in the in- 
tegrand, in this case much the same information is contained 
in the double-integral [F(r p , vi 03 )] as there was in the simpler 
single-integral [f<t>{r,Vd)]. This discovery bodes well for the 
solution of equation M, since, as Pichon & Thiebaut (1998) 
have ably demonstrated, it is quite possible to regularize the 
single-integral case to a point where it is realistically soluble. 
If we are in practice to convert F(r p ,vi os ) into f<j>{r,v$), 
we need the value for the normalization factor, A in equa- 
tion (J^J). On the basis of simple geometric arguments, one 
might expect A to be proportional to r p . However, the exact 
value A cannot be determined in any simple way. Indeed, 
one would expect the value to depend on the DF of the 
disk, the very thing that we are trying to derive. However, 
as we shall see below, the similarity between the shapes of 
F(r p ,vi os ) and f^(r p ,v^) does form the basis for a useful 
iterative scheme for recovering the DF, which is completely 
insensitive to the form adopted for A. 

In order to describe this scheme, we define some simpli- 
fied notation. Let A~ x be the Abel inversion operator of 
equation (n) and V the line-of-sight projection operator in 
equation (|). In this notation, the line-of-sight velocity dis- 
tribution, the DF of the system and the tangential velocity 
distribution are related as follows: 



F(r p ,v los ) =Pf(E,L) =V(A~ 1 U). 



(16) 



We can now present an iterative algorithm which computes 
velocity distributions that converge towards the true tan- 
gential velocity distribution. As discussed above, the part of 
the tangential velocity distribution above the rotation curve 
is sufficient to uniquely specify the DF, and so for the re- 
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mainder of this section we consider only this part of the 
distribution. Let /i n ' be the n th approximation to the tan- 
gential velocity distribution. As discussed in the previous 
paragraph, /i ^ = A _1 F would be a good zeroth-order es- 
timate if we knew the value for A. However, even if A is 
known, this estimate is only an approximation, so A~ x 

will differ somewhat from f(E,L) and hence V(A~ 1 (f^)) 
will have systematic residuals from the observed form of F 
that we started with. As a next iteration, we correct for 
these small residuals by subtracting them from our assumed 
form for the tangential velocity distribution, so that 



fT +1) = fV- A ~ 



nA-\fi n) ))-F 



(17) 



and iterate until convergence is achieved. If this scheme con- 

jM It therefore follows 



verges, we know that 
from equation (|17j) that 



(18) 



Note that, with this formulation, if the iteration converges 
then the result will be totally independent of the adopted 
form for A, so our ignorance as to the exact value for this pa- 
rameter is not a problem - A acts as a relaxation parameter 
in the iteration, which may affect the speed of convergence, 
but not the final answer. Thus, as the iteration converges, 
from equation we have — > fy. One can then ob- 

tain the distribution function by a simple Abel inversion, 
f(E,L) = A- 1 U. 

To implement the Abel inversion operation, A 1 , one can 
draw on the sophisticated techniques developed by Pichon 
& Thiebaut (1998) in their analysis of non-edge-on disks. 
However, this paper is primarily concerned with establishing 
whether the iterative scheme described above converges. For 
this purpose, we can consider model data in which the noise 
level is low, and apply the computationally-simpler approach 
given by equation M. Numerically, at each step, interpola- 
tions of the estimate of the DF, its projected velocity dis- 
tribution and the approximation of the tangential velocity 
distribution /i n ' are computed. We constrain the tangential 
velocity distribution and distribution function estimates to 
be positive. A satisfactory approximation of the DF is found 
when the difference between the LOSVD F and the projec- 
tion of the approximate DF is small everywhere in a large 
region of the (r p ,vi os ) plane. The resolution of the inter- 
polation grid ultimately limits the accuracy with which the 
DF can be recovered with this iterative scheme, eventually 
resulting in noise amplification. 



4 MODEL TESTS OF THE ALGORITHM 

To test the above proposed algorithm, we have constructed 
the line-of-sight velocity distribution as a function of pro- 
jected radius for a moderately-realistic edge-on disk model. 
The DF was taken to be of the form 



/ = .g(L)exp [—/3(L)(E — E C {L))\ , 



(19) 



where g(L) and /3(L) are two arbitrary functions of the an- 
gular momentum and E C (L) is the energy of the circular 
orbit of angular momentum L. Similar forms of DFs were 
described in previous work by Binney (1987) and Kuijken & 
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Figure 1. Top: line-of-sight velocity distribution F(r p , vi oa ) (left) 
and tangential velocity distribution (right) fj,[r, tu); the line 
shows the adopted rotation curve. Bottom: plot of F(r p ,Vi os + 
"c(^p)) and /^(r, v$ + v c (r)) for the region above the rotation 
curve. 



Tremaine (1991). For these tests, we adopt functional forms 
of 



/^)=exp(£) 



and 



g(L) = 







if L > 
otherwise. 



(20) 



(21) 



To complete the model, we must also specify the gravita- 
tional potential. In order to match the flat rotation curves 
found in real galaxies, we adopt the softened isothermal 
sphere potential, 



*(r) 



In 1 + 



(22) 



The upper panels of Fig. [l] show the LOSVD as a function of 
projected radius and the tangential velocity distribution as 
a function of radius for this model. The model does, indeed, 
look very similar to the LOSVDs seen in the stellar kine- 
matics of real edge-on disk galaxies (e.g. Kuijken, Fisher & 
Merrifield 1996). The lower panels show the parts of these 
two functions that lie "above" the rotation curve. For the 
reasons discussed in Section |^, these two plots appear simi- 
lar, allowing us to use the bottom right panel as our initial 
approximation for the bottom left panel. 
We have applied the algorithm described in Section | to a 
model of the form given by Eqs. @, © and ^ 

for a variety of sets of parameters. Figured illustrates the 
case for parameters Lq = 1.6, ro = 0.2, vo = 2.0. As this 
figure shows, the iterative algorithm recovers the DF very 
closely. In fact, convergence occurs in a small number of 
steps. The only remaining discrepancies between the "ob- 
served" LOSVDs and those produced by projecting the re- 
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Figure 2. The distribution function and LOSVD for the simple dynamical model described in the text, as recovered by the iterative 
algorithm Top: the true distribution function expressed as f(^j2[E — E C (L)], L) (left), the distribution function recovered by the algo- 
rithm (middle), and absolute difference between the two (right). Bottom: the true LOSVD (left), the LOSVD recovered by the algorithm 
(middle), and absolute difference between the two (right). 



covered DF occur at very small values of r p , due to numerical 
noise in the integration process. 

Thus far, we have assumed that we know the gravitational 
potential in our reconstruction of the DF. Such an assump- 
tion is reasonable if modelling a disk galaxy containing gas 
from which an emission-line rotation curve can be obtained. 
However, such information would not be available for a 
purely stellar system such as an SO galaxy. Further, we are 
also interested in addressing the more general question of 
whether the gravitational potential is uniquely specified by 
the observable stellar kinematics, or whether one can de- 
rive equally-plausible distribution functions using different 
assumptions about the form of the potential. 

It is apparent from Fig. [j] that there is no obvious way to 
estimate the rotation curve from the observed LOSVDs: the 
combination of asymmetric drift in the stellar kinematics 
and the effects of projection along the line of sight means 
that the local circular speed does not correspond to any sim- 
ple property of the stellar kinematics such as the peak of the 
LOSVD or the mean line-of-sight velocity of the stars. We 
are therefore, in principle, free to choose a different rotation 
curve and hence gravitational potential. 



Figures |3| and |4| show what happens in practice if we do 
so. For these model calculations, we have adopted a gra- 
vitational potential that differs only fairly marginally from 
the true form. The iterative process again converges rapidly 
to a plausible DF, which reproduces the LOSVDs exactly 
for \vi os \ > v c (r p ). However, the LOSVDs that one predicts 
from the derived DF for \vi oa \ < v c (r p ) bear little resem- 
blance to those of the original galaxy model. Thus, it would 
appear that the exact form of the gravitational potential 
is tightly constrained by the observations: using just the 
high-velocity kinematics, one can reconstruct the full DF 
consistent with any given gravitational potential, but the 
low-velocity tails of the LOSVDs will only be correctly re- 
produced if the correct potential is adopted. 



5 CONCLUSIONS 

The ultimate goal of dynamical astronomy is the deriva- 
tion of all that there is to know about a galaxy's dynamics 
from its observable kinematics. We are still clearly a long 
way from attaining this "holy grail," but the analysis of this 
paper does provide some cause for optimism. Specifically, 
we have shown how the distribution function of a relatively 
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Figure 3. Illustration of the effects of assuming a slightly erroneous rotation curve. The left panels show the true properties of the 
model, the middle panels shows the dynamical properties recovered by the iterative algorithm with the wrong adopted potential, and the 
right panels show the difference between the two. The top row shows the distribution functions, the middle panel shows the high-velocity 
part of the LOSVD, and the bottom panels show the full LOSVDs. The parameters of the true potential are (ro = 0.2, vq = 2.0) and 
the parameters used for the reconstruction are (ro = 0.2, vq = 2.1) [see equation (pj[)]. 



simple model galaxy can be estimated directly from its ob- 
served kinematics, using a straightforward iterative scheme. 
Further, the redundancy of information in the kinematics 
means that one can readily rule out models in which the 
wrong gravitational potential has been adopted. 

The success of this iterative scheme seems to derive from 
the fact that the information available from the LOSVD of 
an edge-on disk is very similar to that available for inclined 
disks, for which the inversion to the DF is already well es- 
tablished (Merrifield & Kuijken 1994, Pichon & Thiebaut 
1998). Thus, perhaps rather surprisingly, this problem does 
not appear significantly more ill-conditioned than the sim- 
pler case of the inclined disk, even though an extra integral 
is involved. 

The most tempting practical application for this technique 
is the study of edge-on SO galaxies, since such objects are 
fairly pure stellar disks, which are believed to contain little 



by way of obscuration by dust.0 There are, however, still 
several obstacles to such an analysis. For a start, SO galax- 
ies also contain a central spheroidal bulge component, so 
any attempt to reproduce their kinematics must include a 
suitable dynamical bulge model. In addition, even the disk 
components of these systems are not infinitely thin, so they 
probably also obey a third integral of motion. The simplest 
models would treat the dependence on this third integral as 
a separable function, but there is no reason why real galax- 
ies should follow such simple models. Finally, the data that 



* It is worth mentioning that obscuration will affect the high- 
velocity and low-velocity tails of the LOSVDs in different ways, 
so in a galaxy where emission lines allow us to determine the 
potential unambiguously, a mismatch between the model and ob- 
served low- velocity tails of the LOSVD could be used as a measure 
of extinction. 
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one obtains for real galaxies will be noisy. The derivation of 
LOSVDs from the broadening of spectral lines is a funda- 
mentally noise amplifying process, so even good quality data 
will contain significant noise contribution. Further, in the 
case of the algorithm presented in this paper, we principally 
use only the data from the high velocity side of the LOSVD, 
which will not span a large range in velocities. We will there- 
fore have to obtain rather high dispersion spectra to de- 
termine this function, resulting in a further cost in terms 
of noise. However, the advent of 8-metre class telescopes 
means that the high quality data required for such anal- 
yses should soon be routinely available. Also in its favour, 
the method works directly with the observed LOSVD rather 
than the moments calculated from it. Even the second mo- 
ment of the velocity distribution is very sensitive to noise, so 
a method that avoids calculating such moments is likely to 
be relatively robust. We can also draw on sophisticated Abel 
inversion techniques, such as those developed by Pichon & 
Thiebaut (1998), to minimize any noise amplification. 

Clearly, the next step in this project should be to apply 
these techniques to real data. Preliminary DF reconstruc- 



tions based on LOSVDs obtained from observed spectra of 
edge-on SO galaxies indicate that the algorithm presented 
in this paper can be applied to real data with attainable 
noise levels to yield an estimate of the underlying DF; we 
defer discussion of this development to a subsequent paper 
(Mathieu & Merrifield, 2000). It might therefore reasonably 
be hoped that the application of this approach to real kine- 
matic data will provide a robust method for studying the 
detailed dynamics of disk galaxies. 
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